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We give a relation which links emitter and collector junction voltages, 
V eb and V cb , collector current I e , and the total charge Q b of carriers that 
enter through the basG terminal {electrons in a pnp transistor): 

AQV.b/kT) _ (QVcb/kT) 

I r = const ~ 



This relation is valid for high injection conditions, subject only to minor 
restrictions. The significance for device modeling is discussed. 



I. INTRODUCTION 



A basic concept of charge control theory is that the controlled current 
(collector current) equals controlling charge (base charge) divided by a 
transit time. 1 This paper presents an additional relation which links 
base charge and collector current with junction voltages. The validity 
of this relation is subject only to minor restrictions. When these are 
met the relation holds even under high-level injection conditions. 

Section II presents a derivation of the new charge control relation. 
Equation (15) states the principal result. The discussion in Section III 
points out the significance of this relation for bipolar transistor models. 

II. DERIVATION 

Consider a one-dimensional transistor of pnp polarity. The hole 
current density is given by 

j v = quEp - qDp' (1) 

where the symbols have their customary meaning. We shall assume that 
diflusivity D and mobility y. are related through the Einstein relation: 
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Approximation (a): 

kT 

D = — u. 
Q 

Approximation (bj: It is assumed that electric fields are low enough 
for avalanche multiplication of carriers to be negligible. 

Approximation (b 2 ): The velocity-field relation is idealized by the 
field dependent mobility expression: 

M = 



E 



1+& 

v, 

where Mo = qD /k T is the low-field mobility, considered for convenience 
independent of doping, and where v, is the scattering limited velocity. 
Approximation (b x ) places an upper limit on allowable bias. It is known 
(see for example Ref. 2) that D is underestimated at high fields by 
approximations (a) and (b 2 ) and that approximation (b 2 ) yields too 
gradual a transition from low field velocities to the high field saturated 
velocity. 3,4 Nevertheless, approximations (a) and (b) afford significant 
simplifications in the treatment to follow and are retained for that 
reason. To the extent that our final result, equation (15), is affected 
by them, it must be considered approximate. The errors depend on 
bias and doping profile; they are not expected to exceed a few percent 
for typical situations. The error due to approximation (6 2 ) overempha- 
sizes velocity saturation effects and may be alleviated by choice of 
values of v, larger than the final saturation value in high field regions. In 
high-field regions the current is carried predominately as drift current, 
with a carrier concentration that is nearly constant in such regions, so 
that errors in D are of minor consequence. 

Next we define a quantity a(x) which is the ratio of the hole current 
density at position x to the current density j c leaving the collector 
terminal 

a(x) = *&• (2) 

For direct current conditions, considered here, a approaches unity at 
the collector and is 1/a = (1 + /3)//3 at the emitter. For large common 
emitter current gain /3, a differs negligibly from unity. 
We use 

E = -— $' (3) 
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where yp is the electrostatic potential in units of the Boltzman voltage 
kT/q, and consider equation (1) as a differential equation for p(x). Its 
solution, when p is specified at a point Xi , is 

p(x) = p{x,)e* ~ ~q~D o i x a ^ e (U 

- I*- f a(t) I W | fl *<»-* w dt. (4) 

Equation (4) is valid for any pair of points z, and x. Denote by x E and 
x c the outside edges of emitter and collector transition regions, and 
use equation (4) with .r, = x E and x = x c ■ Multiplication of equation 
(4) by e nz) and use* of 

p(.r) = n t c*> iz) -* {1) (5) 

n(x) = n,e* lT) - p '' (I) (6) 

where <p p and <p n are hole and electron quasi-fermi levels in units of the 
Boltzman voltage, yield 



[ XC a(t)n/ U) dt + — f C a(t) \ *'(«) 



c* l,) dt 



We shall now show that the second term in the denominator is negligible. 
The integrals in the denominator obtain the largest contribution from 
the region near x m where \ft(x) attains its maximum value \f/ m . If in the 
second integral we replace a(t) by its value a m at x„ , and if we neglect 
e^ Ufi) and e* ( ' c> in comparison with e ia — all very reasonable assump- 
tions — we obtain for the second integral in the denominator 



Approximation (c): 



D n n< f' c _,„ , ,„,, , _*to dt _, 2o w D.n. & *. 

v. 



For an assessment of the relative magnitude of the terms in the 
denominator of equation (7), consider that in a region of width w the 
potential yp(x) does not differ markedly from \p m ; such region is con- 
ventionally called the "base" of the transistor. Consider high current 
gain, that is, a « a m » 1. Then the value of the first integral is write* ", 
compared with (2D /v,)n i e' l " n for the second. The quantity 2D /v, has 
units of length and is 5r; 200 A for silicon. This length is small compared 



* Equation (6) is defined for later reference. 
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to base widths of today's most advanced transistors and hence we will 
neglect the second term in the rest of this paper. Conceivably, future 
transistors may have narrow enough bases that the term will have to 
be kept. 

Approximation (d): 

^ f X ° a(t) I Mt) I e* lt) dt « ['" a(«)n/ (,) dt. 

v. J IB J XE 

If in equation (7) we were to let v. — > a>, that is, considered the 
carrier velocity to be strictly proportional to the electric field, then 
approximation (d) would be implemented automatically. Note, how- 
ever, that in making approximation (d) we do not imply v, = °o, nor 
do we neglect essential consequences of the finiteness of v, . A low 
value of v, manifests itself in substantial base widening at high currents, 
that is, in influencing ip(t) in the remaining (first) term in the denomi- 
nator of equation (7). In view of the idealized textbook treatments in 
which the minority canier concentration at the base side of the collector 
depletion region is set equal to zero, rather than to a finite value, the 
following statement of equation (7) may be of interest: For low injection 
(that is, for currents sufficiently low that the base width is independent 
of current) the effect of the finiteness of the scattering limited velocity 
on the dc collector current is equivalent to a base widening of 2D /v, . 

We now make the approximation: 

Approximation (e): 

The value of the electron quasi-fermi 
level in the base is constant. 

A gradient in the electron quasi-fermi level in the region where electrons 
are majority carriers would cause appreciable electron current to flow; 
for transistors of reasonable current gain, such currents are negligible. 
Thus, approximation (e) is very reasonable. We denote this value of the 
electron quasi-fermi level by <p nh and divide numerator and denominator 
of the right side of equation (7) by exp (<p„ 6 )- We define the emitter-base 
and collector-base junction voltages by 

kT 
V tb = y Mx E ) - <p nh ] (8) 

IcT 
V cb = y[«P P (x c ) -**]. (9) 
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These voltages differ from terminal voltages by ohmic drops, primarily 
lateral ohmic drops in the base region. The first integral in the denomi- 
nator of equation (7), after it is divided by exp (<p nb ), contains very 
nearly the total area density of electrons. 

Approximation (f): 

f C nMW"-"^ dt = f' C a(t)n(t) dt. 

•>IB "* Z E 

The integrands outside the base region differ, since there the quasi- 
fermi level is position dependent and does not equal ip nb , but the contri- 
bution to the integral outside of the base region is negligible. By de- 
fining an average value (a)„ of a, 

f " a(t)n(t) dt 

(a) KV = ^- (10) 

/ n(t) dt 

We may write expression (/) as 

n< f° a(t)e n,) -"> dt = -^ g b (11) 

where q b is the total charge, per unit area, of those mobile carriers 
associated with the base terminal, that is, electrons in a pnp transistor. 
Equation (7) with approximations (d) and (f) may be written 

(q a D.nV<aU[ e < ' r ' t/tr) - e " r """] „ 

,. = " (12) 

We now change from current and charge densities to current and charge. 
We chose the sign of the collector current according to the convention 
that an electric current entering the device is positive: 

L = -i,A (13) 

Q b = q b A (14) 

where A is the device area. Note that the sign of Q b is such that an 
electric current flowing into the base tends to increase Q b . This is the 
proper sign for charge control theory. Equation (12) can now be written 

(dV.b/kT) AaV c b/kT) 

Ir = C jr 1 (15) 
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with 

C = frgffigs . (16) 

(a)av 

Equation (15) is the principal result of this paper. Note that Q b de- 
pends on bias, and that the form of the bias dependence is governed 
by the doping profile. However, the relation among the quantities 
I c , V eb , Vch , and Q b in equation (15) is independent of the details of 
the doping profile, provided that assumptions (a) through (f) are valid. 

III. DISCUSSION 

In spite of the simple appearance of equation (15) — indeed because of 
it — it provides a powerful tool for transistor modeling. It may be 
written in the form 

I c = -a 21 [e (QV ' b/kr> - 1] + a 22 [e iQVcb/kT) - 1] (17) 

with 

C 

a 22 = —On = 77- (18) 

which is the form of one of the Ebers-Moll equations. 5 But whereas in 
the Ebers-Moll equation the coefficients a 21 and a n are constant, they 
depend according to equation (18) on bias through the base charge Q b 
(C depends on bias only through (d) av which for $ » 1 is nearly unity 
and varies little with bias). It is this bias dependence of Q b which con- 
tains high-injection effects. Thus, use of equation (15) holds promise for 
transistor modeling of improved accuracy. The major bias dependence 
of the collector current is through the exponentials in the numerator 
of equation (15). These are "ideal" exponentials (unity emission coeffi- 
cients) and involve no approximations. The actual modeling is now 
done on Q b in the denominator. A bipolar transistor model using this 
approach will be presented in a later paper. 
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